Numerical distribution functions for seasonal unit root tests

When working with time series data observed at intervals smaller than a year, it is often necessary to test for the presence of seasonal unit roots. One of the most widely used methods for testing seasonal unit roots is that of HEGY, which provides test statistics with non-standard distributions. This paper describes a generalisation of this method for any periodicity and uses a response surface regressions approach to calculate the critical values and P values of the HEGY statistics whatever the periodicity and sample size of the data. The algorithms are prepared with the Gretl open source econometrics package and some new tables of critical values for daily, hourly and halfhourly data are presented.